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ABSTRACT 


>/  An  on-off  control  provides  a  minimum  time  response  for 
missile  control.  For  application  in  missile  control  systems, 
it  is  wasteful  of  control  effort  (due  to  chatter)  to  use  a 
ideal  relay.  Hence  it  is  necessary  to  modify  the  ideal 
relay  into  a  saturating  linear  control.  The  result  was 
almost  the  same  to  that  of  using  the  ideal  relay. 
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I-  I N THO DOCTIO 8 


Guided  missiles  are  classified  into  four  broad  catego¬ 
ries,  depending  on  launch  and  target  position  characteris¬ 
tics.  These  categories  are  (1)  .  air  -to  -  air  (2).  air  - 
to  -  ground  (3).  surface  -  to  -  air  and  (4).  surface  -to  - 
surface.  Each  category  of  the  above  will  employ  one  or  more 
of  the  following  guidance  schemes;  programmed  command, 
liae-of-sight,  lead-angle,  proportional  navigation  homing 
and  inertial.  The  beam  rider  guidance  is  included  in  the 
line-of-sight  guidance.  &  number  of  missiles  also  use  a 
combination  of  these  methods.  For  example,  the  initial  part 
of  the  missile  trajectory  may  use  programmed  guidance  while 
the  terminal  phase  may  use  beam-rider. 

This  thesis  discusses  the  surf ace-tc-air  missile 
controlled  by  on-off,  thrust  vector,  control.  Consideration 
was  given  to  determine  the  effects  of  the  two-level  relay 
and  the  saturation  linear  control.  In  order  to  verify  the 
results,  it  was  tested  by  using  the  type  of  control  for 
three  different  types  of  missile-target  scenarios; 

(1)  .  LOS  command  against  non-maneuvering  target 

(2)  .  LOS  command  against  maneuvering  target 

(3)  .  Pseudo-LOS  command  against  non-maneuvering  target. 

In  chapter  2,  a  discussion  of  a  line-of-sight  guidance 

was  presented  and  a  practical  example  of  it  was  shown  in 
chapter  3.  The  general  concept  of  on-off  control  was 
described  in  chapter  4.  The  simulation  results  of  the  basic 
connand  to  line-of-sight  against  both  a  non-maneuvering  and 
a  maneuvering  target  were  shown  in  chapters  5  and  that  of 
pseudo-LOS  case  was  in  chapter  6.  Finally,  a  discussion  of 
two-level  relay  and  saturating  linear  control  was  presented 
in  chapter  7.  &  table  of  variables  which  were  used  in  this 
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II.  QVgBYIEW  OF  LINS-OF-SIGHI  GOIDAHCE  CONTROL 

A  LOS  system  can  be  called  a  "3-point”  guidance  system 
since  there  is  one  point  which  defines  the  tracker,  another 
the  target  and  a  third  which  defines  the  position  of  rha 
missile.  The  object  of  the  guidance  system  is  to  constrain 
the  missile  to  lie  as  nearly  as  possible  on  the  line  joining 
the  tracker  and  the  target  called  the  Line  Of  Sight  (LOS). 
The  concept  is  simple  and  can  be  implemented  in  many  ways; 
perhaps  it  is  this  apparent  simplicity  which  explains  why 
many  of  the  guided  weapon  systems  as  yet  designed  are  LOS 
system. 

Consider  a  target  flying  straight  and  at  constant  speed, 
and  a  missile  flying  at  a  different  angle  but  constant 
speed,  having  been  launched  when  the  target  occupies  a  posi¬ 
tion  TO  (see  Figure  2.1). 

After  intervals  of  time  of  1,2,3  etc  seconds  the  LOS  is 
shown  as  OT1,OT2,OT3  etc.  Since  the  mrssile  ideally  always 
lies  on  these  lines  the  flight  path  will  be  a  curved  one, 
for  an  approaching  target,  the  curvature  becomes  increas¬ 
ingly  severe  towards  the  end  of  the  engagement.  We  note  that 
the  tangent  to  the  flight  path  at  any  one  point  defines  the 
instantaneous  direction  of  the  missile  velocity.  It  is  seen 
that  the  missile  velocity  vector  will,  in  general,  not  be 
directed  along  the  LOS;  towards  the  end  of  the  engagement  it 
may  be  at  a  considerable  angle  to  it  [Hef.  1]. 

In  an  actual  situation  the  guidance  signals  transmitted 
to  the  missile  are  the  demanded  lateral  accelerations 
(LATA I)  in  two  axes  at  the  right  angles  to  the  beam.  These 
demands  are  resolved  into  missile  axes  within  the  missile. 
An  error  compensation  term  endeavouring  to  keep  the  error 
off  the  beam  ( tf^)  equal  to  zero. 


Figure  2.1  Hissile  Target  Encounter  with  LOS  Guide 

A  basic  geometry  and  a  simplified  guidance  loc 
shown  in  Figure  2.2  and  Figure  2.3. 

Suppose  that  the  cross  range  error  (CRE)  of  Figi 
can  be  measured  either  directly  or  by  means  of  the  e 
difference  between  OT  and  OM,  together  with  some  knc 
of  missile  range  (Rm),  then 

CRE  =  Rm<  )  (2.  1) 

If  this  error  off  the  beam  is  used  as  an  acceleration 
0#  it  needs  some  damping  so  that  good  response  chara< 
tics  are  obtained.  A  dynamic  equation  of  the  form 

CRE  *  G 1  (CRE)  ♦  G2  (CRE)  (2.2) 

needs  to  be  satisfied,  where  G1  and  G2  are  constants, 
necessity  leads  immediately  to  the  consideration 
filtered  error.  In  the  presence  of  noise  on  the  sigh* 


v  V  ' 
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Figure  2.2  Basic  Geometry 


Figure  2.3  Simplified  Guidance  Loop  of  LOS  Guidance 


1  3 


and  hence  on  the  cress  range  error,  CRE,  such  a  filter 
design  is  not  simple  and  becomes  a  compromise  between 
requirements  for  smoothing  the  noise  and  giving  an  adequate 
response  to  a  demand.  Modern  techniques  allow  filters  to  be 
designed  statistically  if  some  knowledge  of  the  noise  char¬ 
acteristics  is  available  or  can  be  assumed.  Figure  2.3  shows 
the  position  of  such  a  filter  ?(s)  in  the  guidance  loop.  It 
includes  a  gain  G,  and  the  acceleration  demand  is 

0  =  F(s)  Rm  (  <5t  “  6m  )  (2.3) 

The  missile  transfer  function  is  represented  by  A(s)  and 
when  the  achieved  acceleration  is  doubly  integrated  and 
divided  by  Rm  it  represents  a  new  measure  of  the  missile 
beam  angle  (tfm)  ,  thus  closing  the  loop  when  differenced  with 
the  target  beam  angle  (dt)  • 

While  this  concept  is  simply  a  LOS  or  beam  riding 
guidance  situation  it  is  by  no  means  as  clear  in  homing  how 
a  guidance  law  can  be  devised  in  the  absence  of  information 
on  missile  and  target  positions  [Ref.  2]. 


Ill-  TYPICAL  BUG AGBME8T  SEQOBHCE 


In  order  to  provide  a  "vehicle"  through  which  to  better 
understand  the  basic  aspects  of  command  to  line-of-sight 
guidance  methodology,  the  engagement  sequence  of  a  short- 
range,  air-defense,  aissile  system  is  described.  The  Boland 
system  was  selected  because  the  general  operational  aspects 
of  the  system  are  available  at  the  unclassified  level 
[Ref.  4]. 

The  entry  of  one  or  more  aerial  targets  into  the  range 
of  the  search  radar  is  indicated  to  the  Roland  vehicle 
commander  by  an  audible  tone.  At  tha  same  time,  a  synthetic 
display  of  the  targets  appears  on  a  screen  to  give  the 
comm aander  all  the  information  needed  to  select  the  most 
threatening  target.  The  screen  images  are  different  for 
friendly  and  enemy  targets.  Also,  the  entry  of  the  target 
into  the  aissile  envelope,  utilizing  target  advanced-range 
computations,  is  indicated  by  a  change  in  the  display.  With 
the  search  antenna  raised  and  the  search  radar  activated, 
target  acquisition  is  possible  even  when  the  vehicle  is  in 
motion. 

There  are  three  modes  of  identification,  friend  or  foe 
(IFF)  interrogation:  automatic,  manual,  and  automatic 

within  a  given  range. 

When  the  commander  has  recognized  a  target  as  hostile 
and  decided  to  engage  it,  he  places  a  cursor  over  the  screen 
image.  This  automatically  brings  the  turret  to  bear  and 
tracking  can  commence  in  either  the  "radar "  or  "optical" 
modes. 

In  the  "radar"  mode,  the  tracking  radar  automatically 
accepts  target  designation  from  the  search  radar,  searches 
for,  locks  onto,  and  tracks  the  target. 


In  the  "optical"  mode,  the  aimer  searches  for  the  target 
in  elevation  with  an  optical  sight.  To  aid  him  an  electronic 
instrument  displays  the  maximum  theoretical  elevation  for 
the  search.  3hen  the  aimer  has  acquired  the  target  in  his 
cross-hairs,  he  keeps  the  target  in  his  sight  by  manipu¬ 
lating  a  control  stick.  This  control  keeps  the  target  prop¬ 
erly  positioned  by  moving  the  turret  in  azimuth  and 
swivelling  a  mirror  in  elevation. 

As  soon  as  the  commander  confirms  that  the  target  is 
within  missile  range,  he  initiates  the  firing  sequence  in 
the  "radar"  mode,  or  authorizes  "optical"  mode  firing 
through  a  command  displayed  in  the  aimer's  sight.  The  aimer, 
then,  can  initiate  the  firing  sequence. 

The  missile  is  guided  by  a  command  to  line-of-sight 
technique.  This  means  that  the  target  is  tracked  optically 
or  by  radar  and  the  deviation  of  the  missile  from  this  line 
of  sight  is  determined  and  corrected  by  a  guidance  command. 
The  commander  may  switch  from  "radar"  to  "optical"  and  back 
again,  as  desired,  even  after  the  missile  has  been  launched. 

Target  tracking  and  determination  of  the  missile's  devi¬ 
ation  from  the  line  of  sight  are  different  for  each  mode.  In 
the  "radar"  mode,  the  guidance  radar  has  two  receiving 
channels.  One  is  used  for  target  tracking  and  the  other  is 
used  to  locate  the  missile  in  the  radar  lobe  through  recep¬ 
tion  of  the  missile's  radio  frequency  beacons.  By  comparing 
these  angles,  an  error  between  the  missile  and  the  target 
line  of  sight  can  be  determined.  In  the  "optical"  mode,  a 
biaxially-stabilized  mirror  is  manually  controlled  to  keep 
the  target  vertically  in  the  aimer's  sight  and  the  turret  is 
rotated  to  the  azimuth  of  the  target  line  of  sight.  An 
infrared  goniometer  is  mounted  to  provide  misile  angle  from 
the  tracker  by  following  flares  mounted  on  the  rear  of  the 
missile.  Then,  a  deviation  of  the  missile  angle  from  the 
target  line  of  sight  can  be  determined. 


Two  groups  of  signals  are  introduced  into  the  command 
computer:  the  velocity  of  the  line  of  sight  in  azimuth  and 
elevation,  and  the  deviation  of  the  missile  from  the  line  of 
sight  in  azimuth  and  elevation.  Based  upon  data  from  the 
liae-of- sight  movement  and  the  angular  deviations  of  the 
missile,  the  necessary  guidance  signals  are  calculated. 

The  guidance  signals  are  relayed  to  the  missile  by  a 
transmitter  with  highly  directional  characteristics.  The 
command-transmitting  antenna  is  slaved  to  the  missile  angle 
in  both  azimuth  and  elevation.  It,  therefore,  is  trained  on 
the  missile  continuously. 

The  side  forces  required  for  missile  course  corrections 
are  produced  through  deflection  of  the  exhaust  jet  of  the 
sustainer  motor  by  spoilers  at  the  rear  of  the  missile 
(thrust- vector  control)  . 

When  the  missile  reaches  the  point  of  impact  with  the 
target,  the  warhead  is  detonated  by  either  percussion, 
contact  fuse  or  the  radio- frequency ,  proximity  fuse.  The 
warhead  consists  of 'a  radial-effect,  multiple-fragmentation 
charge. 

Figure  3.1  presents  an  operational  schematic  of  the 
basic  Roland  missile  system  operation. 

The  computer  simulations  contained  herein  are  generic  in 
nature  within  the  command  to  line-o f-sight  guided-missile 
type  and  have  only  reasonable  estimates  of  missile  capabil¬ 
ities  introduced.  This  ensures  unclassified  results.  At  the 
same  time,  the  simulations  are  of  sufficient  complexity  to 
properly  weigh  the  relative  merits  of  the  guidance 
variations  discussed  [Ref.  31. 


17.  OB-QFP  (BANG-BANS)  COHTBOL 

As  discussed  before,  LOS  guidance  maintains  a  missile 
position  on  the  LOS.  Usually  missile  position  has  a  cross 
range  error  (CRE)  and  we  want  to  raduce  this  error  to  zero 
in  the  ainiaua  tiae.  This  kind  of  problem  can  be  solved  by 
using  the  on-off  control.  The  basic  concept  of  this  is  that; 

Given  a  system  for  which  the  drive  is  limited  (has  a 
maximum  or  saturation  value) ,  the  fastest  response  is 
obtained  if  maximum  forward  drive  is  applied  at  t  *  0, 
and  is  reversed  at  a  proper  instant  t  =  1 1  so  that 
deceleration  under  maximum  reverse  drive  reduces  the 
velocity  to  zero  at  precisely  tha  command  value  of  the 
output.  The  drive  is  then  set  to  zero. 


Tha  ideal  relay  permits  only  two  conditions;  full 
acceleration  and  full  deceleration  [Ref.  5]. 

Proa  the  Bang-Bang  control  law,  we  can  derive  tha 
switching  function  which  makes  the  arror  go  to  zero  by  using 
the  proper  switching  time.  Prom  Newton's  second  law; 
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at  t  * 

0  ,  ere 

a  o  and 

kl  ■ 

0.  Therefore 

CRE  « 

d 

— —  (CRE) 
dt 

■  U  t 

(4.  1) 

CRE  ■ 

CRE  dt 

a  g  t2 

♦  k2 

(4.2) 

19 

v. 


From  the  equation  4.1 


t  a  CRE  /  U 

t2  =  (CRE  /  U)2  (4.3) 

Substitute  equation  4.3  into  equation  4.2 

CRE  a  (0  /  2)  (CRE  /  U)2  ♦  k2 

=  (CREJ2  /  (2  0)  ♦  k2  (4.4) 

where  lc2  is  integration  constant.  But  if  we  apply  a  full 
deceleration  at  the  halfway  point,  the  equation  4.4  becomes 

F  a  (CRE  |  CR E|  )  /  (2  0)  ♦  CRE  (4.5) 

and  is  called  the  ERROR  FUNCTION.  U  will  be 

U  a  +  G 

or 

U  =  -(G)  SIGN(F)  (4.6) 

Equations  4.5  and  4.6  represent  the  SNITCHING  FUNCTION  which 
makes  the  error  go  to  zero  in  the  minimum  time.  The 
switching  function  and  the  block  diagram  of  the  on-off 
controller  are  depicted  on  the  Figures  4.1  and  4.2.  And  we 
can  obtain  the  cross  range  error,  CRE,  by  doubly  integrating 
U  with  the  initial  condition  of  CRE.  we  have 

CRE  »  U  dt  ♦  CRE (0) 

CRE  »  CRE  dt  ♦  CRE  (0)  (4.7) 

The  simulation  results  of  these  equations  are  giver,  on 
Figures  4.3  through  4.7  and  the  computer  program  is  attached 
(see  Appendix  B)  . 


Figure  4.1  Parabolic  Switching  Fanction 


r 
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Figure  4.2  Bloch  Diagram  of  ON-OFF  Controller 


Figure  4.5  F  versus  Time 


CRE 


CRED  VS.  CRE 
SWITCHING  FUNCTION 


V.  BA- I C  COMMAHD  TO  LINE-OF- SIGHT  SIMULATION 

A.  SCENARIO 

The  engagement  was  designed  with  the  ground  tracker  and 
missile  launching  unit  located  at  taa  origin. 

The  target  was  flown  accross  the  first  quadrant  from  a 
position  4000  meters  on  the  x-axis  and  1000  meters  on  the 
y-axis  (4000,1000).  The  velocity  vector  of  the  target  was 
parallel  to  the  x-axis  and  magnitude  was  250  meters  per 
second. 

Since  most  missiles  need  a  few  seconds  of  boost,  the 
missile  is  not  controlled  during  this  time.  He  assumed  that 
the  missile  was  controlled  after  one  second  from  the  firing 
time  and  controlled  by  PROGRAMMED  GUIDANCE  up  to  this  time. 
After  the  time  of  missile  ••capture" ,  the  missile  was 
controlled  by  the  on-cff,  TVC  method  with  the  LOS  guidance 
law.  The  simplified  flow  chart  is  shown  in  the  Figure  5.1. 

In  order  to  simplify  the  problem,  we  assumed  that: 

1)  the  velocity  vector  of  missile  ,  Vm,  was  parallel  to 
the  LOS  between  the  target  and  origin  and  the  magni¬ 
tude  of  7m  was  constant,  503  meters  per  second; 

2)  the  LATAX  was  applied  to  the  missile  at  right  angles 
to  the  LOS.  This  was  a  reasonable  assumption  for 
this  kind  of  missile.  So  the  angle  <P  +  in  the 
Figure  2.2  is  almost  same  to  angle 

3)  the  measurement  noise  was  zero  so  we  omitted  the 
filter ,  F(s)  ; 

4)  the  magnitude  of  LATAX  was  150  meters/second  which 
was  about  15  Gs. 

The  geometry  depicted  in  Figure  5.2  summarizes  the  geometric 
situation. 


25 


Figure  5. 


1  simplified  Flow  Chart  of  Basic  LOS  Coea 
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Figura  5.2  Geometry  of  Basic  LOS  Guidance 

For  mathmatical  convenience  of  simulation,  we  need  to 
define  the  sign  of  the  CRE  and  the  LATAX  as  follow; 

♦  (CRE|  :  When  the  missile  position  is  upper-side  of  LOS 

-  |CRE |  :  when  the  missile  position  is  lower-side  of  LOS 

♦  1 0|  :  when  the  LATAX  is  upward  direction 

-  I Q |  ;  when  the  LATAX  is  downward  direction 

This  sign  was  based  on  the  positive  which  is  defined 

when  is  greater  than  6^  • 

B.  PRO GRAB BED  G 01 DANCE  PHASE 

Since  the  major  emphasis  of  this  paper  was  on-off 
control,  we  assumed  that  the  missile  flew  along  the  LOS 
during  the  programmed  guidance  phase.  3ut,  in  a  practical 
situation,  there  is  some  cross  range  error  which  is  occured 
by  disturbances  such  as  wind,  propulsion  system  and 


autopilot  time  delay, etc.  Hence  we  made  ini tializatior. 
errors,  and  the  on-off  control  started  with  these  errors. 

C.  OH-OFF,  THRUST  7ECT0B,  MISSILE  CONTROL 

The  detail  of  the  on-off  control  was  discussed  before, 
hence  we  applied  this  to  the  LOS  guidance  scheme.  The  block 
diagram  of  this  system  is  depicted  in  the  Figure  5.3 

(Ref.  3]. 

In  order  to  determine  the  CRE,  the  tracker  estimates  the 
missile's  range  (Rm)  ,  by  the  elapsed  time  of  flight  and  the 
missile's  velocity  profile.  The  program  of  this  simulation 
is  attached  in  Appendix  C. 


Missile 

Position 

Geometry 


TRACKER 


>mtJ  Control  1  j  F( 
1  Function  fl 


MISSILE 


Figure  5.3  Block  Diagram  of  the  Basic  LOS  Command 


D.  SIMULATION  RESULTS 

Figure  5.4  shows  the  missile  and  target  geometry  in  X-Y 
plane.  The  missile  intercepted  the  target  at  the  point 
A(2605,1000)  with  the  almost  zero  miss  distance. 


Figure  5.5(a)  shows  the  distance  between  target  and 
aissile  versus  tine.  The  distance  decreased  linearly  and 
neared  zero  at  the  time  at  5.58  seconds. 
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Pigure  5.4  The  Basic  LOS  Coaeand 

Figure  5.5(b)  shows  the  CRB  versus  time.  The  initial 
cross  range  error  at  the  missile  capture  time  one  second 
after  firing  was  about  50  meters.  Since  the  CEE  at  the 
"missile  capture  time"  was  positive,  the  CRE  increased 
initially.  So  the  maximum  CRE  was  about  58.2  meters  at  the 
time  1.330  seconds.  Then  it  decreased  to  almost  zero  meter 
at  2.55  seconds.  To  get  a  faster  response,  we  should 
increase  the  magnitude  of  the  LAI&X.  Me  should  note  here 


that  the  CRE  does  not  maintain  zero  value  because  we  did  not 
consider  the  target  motion  terms  in  this  phase  of  the  simu¬ 
lation  (Ref.  2].  So  the  missile  had  some  small  cross  range 
error  and  the  BANG- BANG  controller  had  tried  to  reduce  this 
error  in  a  chatter-mode. 

Figure  5.6(a)  shows  the  CRE  versus  time.  Figure  5.6(b) 
shows  the  CRE  versus  CRE.  As  we  expected  this  curve  followed 
the  SNITCHING  FUNCTION  as  shown  in  figure  4.4.  Figures 
5.7(a)  and  5.7(b)  show  the  F  versus  time  and  the  0  versus 
time. 

This  program  was  tested  using  maneuvering  targets  and 
the  results  were  almost  same  except  the  impact  position.  The 
results  cf  this  simulation  were  shown  on  the  Figures  5.8 
through  5.11  and  the  program  is  attached  in  Appendix  D.  The 
comparison  of  these  simulations  is  summarized  in  Table  I. 


Pigure  5.8  The  Maneuvering  Target 


CRED  VS.  TIME 
WITH  MANEUVERING  TARGET 


The  Haneuvering  Target 


CRED  VS.  CRE  ERROR  CFJ  VS.  TIME 

WITH  MANEUVERING  TARGET  g  WITH  MANEUVERING  TARGET 


The  Maneuvering  Target 


Figure  5.11 


versus  Time  for  a  maneuvering  Target 


TABLE 

I 

The  Basic  LOS  command 

Simulation  Result 

NON-MANEUVRING  TARGET 

MANEUVERING  TARGET 

.  time  (control)  *  1.0  sec 

CRE  (0 )  »  49.910(m) 

CRED(O)  =  49. 8 32  (m/sec) 

1.0  sec 

49.  907  (m) 

49.891  (m/sec) 

.time (max. CRE)  1.33  sec 

CRE  (max)  *  58.  184(m) 

1 . 3  3  (sec) 

58.201  (m) 

. time (intercept) =  5.58  sec 
(Xm  ,  im)  =  (2504.7,999.  88) 
CRE  *  4.3344  E-5 

distance®  0.348  94 (m) 

5.61  sec 
(2597.  9,  1057.7) 
-6.5448  E-6 

0.47889  (m) 

71.  PS SDDO-LOS  COMMAND  SIMDLATIOH 


The  guidance  scheme  cf  the  lead  angle  command  is  almost 
the  same  as  that  of  the  basic  LOS  command.  Instead  of  the 
tracker- to- tar  get  line-of- sight  this  guidance  scheme  uses 
the  tracker-to-estimated  impact  point  and  is  called 
’'synthetic  line- of- sight"  (SLOS)  ,  or  "pseudo  line-of-sight". 
The  missile  is  controlled  to  fly  along  this  pseudo  line-of- 
sight.  The  block  diagram  of  this  system  is  easily  modified 
from  that  of  the  basic  LOS  and  is  shown  in  Figure  6.  1. 

The  estimated  impact  point  at  the  instantious  time  is 
calculated  by  using  the  "time  to  go"  (Tg)  and  the  "closing 
velocity"  (Vc|  between  the  target  and  the  missile. 
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Position1 
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Velocity 


T.V.C. 


SuiSBlpil 


Missile 

Position 

Geometry 
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The  "closing  velocity"  and  "time  to  go"  are  calculated  as 
follow: 

Vc  =  ((Vtx-Vmx)1  ♦  (Vty-Vmy*)^ 

Tg  =  (distance  between  target  and  missile)  /  Vc 
=  f(Xt-Xm)2  +  (Yt-Ym)2)  *  /  Vc 

The  missile  goes  to  the  impact  point  directly.  The 
simulation  result  is  almost  same  as  in  the  basic  LOS  case. 
On  other  hand,  this  guidance  scheme  is  poor  in  a  ECU  situ¬ 
ation.  In  order  to  compare  the  results  we  used  the  same 
data  as  that  of  the  basic  LOS  command.  These  are  shown  in 
Figures  6.2  through  6.5  and  the  summarized  results  are  shown 
in  Table  II.  The  computer  program  is  attached  in  Appendix  E. 


TABLE  II 

The  Pseudo  LOS  Command  Simulation  Result 


,  time  (control) 
CRE(Q) 


.time  (MAX.CRE) 
CRE  (max) 

.time  (intercept) 
(Xm,  Ym) 

*B,rt 

miss-distance 


1.0  sec 
^9.912  (m) 

49.829  (m/sec) 

1.33  sec 
58.  184  (m) 

5.58  sec 
(2604.7,  999.92) 
(2605.0,  1000.0 
2.  9572  E-5(m) 

0.  35137  (m) 


CRED  VS.  CRE  ERROR  (F)  VS.  TIME 

PSUEDO  -  LOS  COMMAND  g  PSUEO0  -  LOS  COMMAND 


Figure  6.4  The  Pseudo-LOS  Coeeand 


71 1.  SIMULA  TIP  8  S  THO-LBVBL  RELAY  AND  SATURATION  CONTROL 

The  LOS  guidance  with  an  "ideal"  relay  has  been 
discussed.  In  this  chapter,  the  effect  of  the  different 
types  of  lateral  acceleration  demand  are  discussed.  In 
order  to  compare  the  results  with  the  previous  simulations, 
the  same  parameters  were  used. 


A.  THO-LEVBL  RELAY 

The  large  magnitude  of  the  LATAX  makes  a  fast  response. 
But  in  the  case  of  small  CRE,  a  smaller  magnitude  of  LATAX 
is  needed.  This  idea  was  developed  in  a  "two- level"  relay  as 
shown  on  the  Pigures  7.1(a;b).  The  shaded  area  on  Figure 
7.2(b)  shows  the  region  of  a  lower  level  of  LATAX  in  the 
"CRE  verse  CRE"  phase  plain.  It  provided  the  minimum  over¬ 
correction.  The  computer  programs  were  easily  modified  by 
adding  one  statement, 

IF  ((|CRE|  ♦ICRSDI)  .LT.  M)  G  =  3/(N1/N2). 

He  used  the  values  150  m/seca  for  N1  and  15  m/seca  for 
N2  and  1.0  for  M  in  the  simulations  of  the  basic  LOS  command 
and  the  pseudo  LOS  command.  The  results  were  almost  the  same 
as  the  previous,  except  in  the  figure  for  "U  versus  time". 
Table  III  summarized  these  simulation  results  and  Figures 
7.2,  7.3  and  7.4  show  the  "U  versus  time"  of  each  case.  The 
programs  were  attached  in  Appendix  P,  G  and  H. 
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Two-Level  Belay  Control  Result 
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ft 

UNIT 

NON-HVR 

aVR-TGT 

PS  3  UDQ-LOS 

^■>: 

.  time  (ccntro  1) 

CHE  (0) 

CRED(O) 

1.0 

49. 910 

1 . 0 

49.907 

1.0 

49.910 

-  %. 

49.  832 

49.891 

49.828 

.  time  (MM.  CHE) 

1.33 

1.33 

1.33 

r.  -•/• 

r  V#  . 

CHE  (max) 

58.  184 

58.201 

58. 184 

,\\\\ 

.  time  (intercept) 

5.58 

5.61 

5.58 

r.  •' 

Xm 

2604.7 

2597.9 

2604.7 

7m 

999.88 

1057.7 

999.92 

Xt 

2605.0 

2597.5 

2605.0 

Yt 

1000.0 

1057.5 

1000.0 

r.-  — 

CHE 

-2. 63E-6 

-3.63E-6 

5.03E-8 

.  ** 

miss-distance 

0.  34894 

0.47889 

0.35137 
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Figure  7.2  0  versus  rise  for  the  Basic  LOS  Guidance 

B.  SATURATING  LINEAR  CONTBOL 

In  the  previous  section  the  two-level  relay  was 
discussed.  The  "sautrating  linear  control"  as  depicted  in 
the  Figure  7.5(a)  and  7.5(b)  was  also  studied.  The  shaded 
area  on  the  Figure  7.5(b)  shows  the  region  of  linear  control 
in  the  "CBE  versus  CRE"  phase  plane.  The  computer  programs 
were  easily  aodified  by  adding  one  statement, 

IF  (ABS  (F)  .LE.N)  U  =  -3  *  F  /  N. 

The  value  of  "H"  determines  the  linear  region  for  F. 

« 

The  Figures  7.6  through  7.8  show  the  simulation  results  of 
the  basic  LOS  command  against  the  non-maneuvering  target 
case  for  "H"  equal  1,  5  and  10.  Hh9n  choosing  the  value  "N" 
equal  to  "one",  the  intercept  time  and  miss  distance  are 
almost  the  same  as  the  counterpart  of  the  ideal  relay  case. 
Hence  the  saturating  linear  control  can  be  used  in  practice 
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Figure  7.3  0  verse  Tiae  against  the  HVB  Target  with 

T vo-Level  Belay 


TIME 


instead  of  the  ideal  relay  by  choosing  a  proper  value  of 
"ft".  The  summarized  results  are  in  the  Table  IV. 

Figures  7.9  through  7.11  show  the  results  of  the  maneuvering 
target  case  and  Figures  7.12  through  7.14  show  the  results 
of  the  pseudo-LOS  case.  These  programs  are  given  in 
Appendix  I,  J  and  K . 


Figure  7.5  Linear  Switching  Belay 


43 


T  ABLE 

I? 

Saturating 

Linear  Control  Result 

(H=1) 

QUIT 

NON-MVR 

MVR-TGT 

PS  SO  DC-LOS 

.time  (control) 
CRE(O) 

CRED  (0) 

1.0 

49. 910 
49.  832 

1 . 0 

49.907 

49.891 

1.0 

49.912 

49.828 

.time  (MAX.  CRE) 

CRE  (max) 

1.33 

58.  184 

1.33 

58.201 

1.33 

58. 184 

.  time  (intercept) 

Xm 

Ym 

Xt 

Yt 

CRE 

miss-distance 

i 

5.58 

2604.7 

999.84 

2605.0 

1000.0 

-.0349 

0.  3  507 

5.61 

2598.0 

1057.6 

2597.5 

1057.5 
-.0702 
0.4841 

5.58 

2604.7 

999.88 

2605.0 

1000.0 

-.0349 

0.3492 

4 

4  9 


The  Basic  LOS  Coaaand  (Saturation) 


Figure  7.7  The  Basic  LOS  Coamand  (Saturation) 
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Pigore  7.8  The  Basic  LOS  Comaanfl  (Saturation) 


Figure  7.9  Haneuwering  Target  (Saturation) 


TIME 


Maneuvering  Target  (Saturation) 


Figure  7.12  Pseudo-LOS  Coauand  (Saturation) 


7111.  COHCLOSION 


The  comparision  of  results  for  the  ideal  relay,  two- 
level  relay  and  saturating  linear  control  against  the  basic 
LOS  command  and  pseudo-LOS  command  against  a  non-maneuvering 
target  and  maneuvering  target  are  provided  in  the  Tables  7, 
VI  and  VII.  These  simulation  results  clearly  demonstrate 
that  "on-off"  control  of  a  missile  is  highly  desireable  and 
that  "saturating  linear  control"  of  a  missile  has  little 
adverse  effects  compared  to  an  "ideal  relay"  control. 


TABLE 

7 

Comparison  of 

the  Basic 

LOS  Command 

Results 

UNIT 

IDEAL 
REL  AY 

TWO- LEVEL 
RELAY 

SATURATION 

CONTROL 

.  time  (control) 

CRE  (0) 

CRED(O) 

1.0 

49. 910 
49. 832 

1.0 

49.910 

49.832 

1.0 

49. 910 

49.832 

.time  (MAX.  CRE) 

CRE  (max) 

1.33 

58.  184 

1.33 

58.184 

1.33 

58. 184 

.  time  (intercept) 

Xm 

Ym 

CR  E 

miss-distance 

5.58 

2604.7 

999.88 

4.33E-5 

0.34894 

5.58 

2604.7 

999.88 

4.33E-5 

3.34894 

5.58 

2604.7 

999.88 

4.33E-5 

0.34894 

_ i 

TABLE  71 

Coaparisoa  of  the  Maneuvering  Target  Results 


UNIT 


.time  (control) 
CRE  (0) 

CRED(O) 

.  time  (MAX.  CRE) 
CRE  (max) 

.time  (intercept) 
Xm 
Tm 
CR  E 

miss-distance 


IDEAL 

THO- LEVEL 

SATURATION 

RELAY 

RELAY 

CONTROL 

1.0 

1 . 0 

1.0 

49. 907 

49.907 

49.907 

49. 891 

49.891 

49.  891 

1.33 

1 . 33 

1.33 

58.  201 

58.201 

58. 201 

5.6  1 

5.61 

5.61 

2597.5 

2597.5 

2597.5 

1057.5 

1057. 5 

1057.5 

-6. 54E-6 

3 . 63E-6 

-0.0702 

0.47889 

0.47889 

0. 4841 

TABLE  711 


Comparison  of 

the  Pseudo- 

L35  Command 

Results 

UNIT 

IDEAL 

RELAY 

TWO- LEVEL 
RELAY 

SATURATION 

CONTROL 

.  time  (control) 

CRE  (U) 

CRED(O) 

1.0 

49. 912 

49. 828 

1.0 

49.910 

49.828 

1.0 

49.912 

49. 828 

.time  (MAX.  CRE) 

CRE  (max) 

1.33 

58. 184 

1.33 

58.184 

1.33 

58. 184 

.time  (intercept) 

Xm 

7m 

CH  £ 

miss-distanca 

5,58 

2604.7 

999.92 

2.96E-5 

0.  35137 

5.58 

2604.7 

999.92 

5.03E-8 

0.35137 

5.58 

2604.7 

999.92 

2.96E-5 

0.35137 

&PPBSDII  A 
VARIABLES  LIST 


DIAGRAH 

VARIABLES 

COHPOTER 

VARIABLES 

NOON  DESCRIPTION 

CRE 

CRE 

cross -range- err or 

CRE 

CRED 

rate  of  cross-range- error  change 

P 

P 

error  function 

G 

G 

magnitude  of  lateral  acceleration 

0 

0 

missile's  lateral  acceleration 

<$t 

SIGT 

angle  between  the  LOS  to  target 

and  X-axis 

<5m 

SIGH 

angle  between  the  beam  to  missile 

and  X-axis 

</nt 

SIGMT 

angle  difference  between  a  and 

Vc 

VC 

crossing  velocity 

tg 

TG 

tiae  to  go 

t  (control) 

TCON 

beginning  tiae  of  on-off  control 

SLOS 

SLOS 

synthetic  line-of- sight 

Xm 

XH 

X-coordinate  of  missile  position 

Ya 

YH 

Y-coordiaate  of  missile  position 

It 

XT 

X-coordinate  of  target  position 

It 

YT 

Y-coordinate  of  target  position 

Va 

VH 

velocity  of  the  missile 

vt 

VT 

velocity  of  the  target 

Vax 

VH X 

X-coaponant  of  missile's  velocity 

Vay 

VHY 

Y-coaponant  of  missile's  velocity 

Vtx 

VTX 

X-compoaant  of  target's  velocity 

vty 

VTY 

I-coaponant  of  target's  velocity 

AP PEN  PI I  B 

PROGRAM  OF  THE  SWICHING  PUNCTION 


TITLE 

TITLE 

TITLE 

INTGER 

CONST 

INITIAL 


BANG-BANG  CONTROL 
SWITCHING  FUNCTION 
*  YE  UN, 

NPLOT 
N PLOT® 1 


J.Y, 


CRE  =  1.0 
CRED  *  0. 
ACRED  =  0. 

********** 

DERIVATIVE 

********** 

NOSORT 


G  =  1.0 
F  =  CRE 


SAMPLE 


_  ♦  <CRED*ACRED) /(2*G) 

U  =  -G  *  SIGN ( 1 . ,P) 

CRED  ®  INTGRL(0.,U) 

CRE  =  INTGBL  (CRE,  CRED) 

ACRED  ®  ABS  (CRED) 

IF  (CRE.  LE .  0.0)  CALL  ENDJOB 


CALL 

CALL 

CALL 

CALL 

CALL 


DRWG 
DR  WG 
DRWG 
DRWG 
DRWG 


1  ,  1, TIME,  CRE) 

2.1,  CRE,  CRED) 

3.1,  TIME,  Ff 

4.1,  TIME, CRED) 

5.1,  TIME, U) 


TERMINAL 

CALL  ENDRW  (NPLOT) 

CONTRL  FINTIM®  2.  1,DELT=0.01  ,DELS=0.01 

PRINT  - - 

END 
STOP 


0.  1  ,G,F,U, CRE, CRED 
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APPENDIX  C 

PROG BAN  OP  THE  BASIC  COS HAND  TO  LOS 


TITLE  BASIC  COMMAND  TO  L.O.S 

TITLE  WITH  IDEAL  RELAY 

TITLE  *****  YEUN,  J.Y.  ***** 

INTEG  RKSFX 

INTGER  N P LOT  .KILL 

CONST  N PLOT= 1 .TCO  Na1 . 0 

CONST  VM=530. ,VT= 250, ?I= 3. 1 41 593,  KILL=0 

INITIAL 

XTO  =  4300. 

YTO  =  1  000. 

GANT  =  PI 
F  =0. 

* 

***************************  ********************************* 
DERIVATIVE 

************************************************************ 

NOSORT 

************************* 

***  TARGET  PARAMETERS  *** 

************************* 

VTX  =  VT*COS(GAMT) 

VTY  »  VT*SINjGAMTi 
XT  =  VTX*TIME  ♦  XTO 
YT  =  VT Y*TI  ME  ♦  YTO 
SIGT  =  ATAN2(YT,XT) 

IF  (II  ME.  GE. TCO N)  GO  TO  50 
**** ********** *********** 

***  MISSLE  PARAMETERS  *** 

************************* 

*  PROGRAMMED  GUIDANCE  ************************************** 

* 

SIGM  a  SIGT  +0.  1 
VXM  =  V  M  *  COS  (SIGM) 

VYM  =  VM  *  SIN  (SIGM 
XM  a  I  NTGRL  (0.  .  VX  M 
YM  a  INTGRL(0.  .VYM 
RM  *  SQBT(XM**2  ♦  YM**2) 

SIGHT  a  sidH  -  SIGT 
CRE  =  RM  *  SIN  (SIGHT) 

CRED  a  DERI  V  (0  .  .CR  E) 

GO  TO  200 

***************************** 

*  ON-OFF  GUIDANCE  (BANG-BANG  CONTROL)  ********************* 
***************************** 

* 

50  CONTINUE 


G  *  150. 

ACRED  «  ABS  (CRED) 

F  »  CRE  ♦  (CRED*ACRED)/(2*3) 
0  *  -G  *  S IGN ( 1 .  .  F ) 

CRED  *  INTGRLlO.O,  U) 

CRE  «  INTGRL(CRE.CRED) 

RM  *  V  H*T  I  ME 
A3  »  CRE/RM 
SIGHT  *  ARS IN ( A3) 

SIGM  »  SIGT  ♦  SIGHT 
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x n  =  aa  *  cos (sigm) 

YM  =  R  M  *  SIN  (SIGH) 

***************  ************  ******************  ************** 
200  CONTINUE 

* 

***************************  *************************  ******* 
***  MISSION  RESULT  ***  KILL  =  0  ;  TGT  MISSED 

**********************  KILL  =  1  ;  TGT  DESTROYED 

XDIST  =  XT-XM 
YDIST  =  YT-YM 

DIST  *  S  QRT  (XDIST  **2  ♦  YDIST**2) 

IF  (DIST  .  LE.5)  KILL  =  1 

IF  (DIST  .GT.51  KILL  =  0 

IF  (XM  .GT.  (XT +30  )  )  C ALL  ENDJOB 

* 

SORT 

*************************************************  ********** 

********  OUTPUT  AND  PLOT  CONTROL  CARD  it******************** 
********************************************************* ** 


SAMPLE 


TERMINAL 


CALL  D RWG  ( 
CALL  DR WG  I 
CALL  D RWG l 
CALL  DRUG  i 
CALL  D  RWG 
CALL  DRWG? 
CALL  DRWG i 
CALL  DRWG  < 


CALL  ENDRW  (NPLOT) 

CONTRL  FINTIM=6 . 0, DELT=0.00 1,DELS=0. 003 
PRINT  0.  005,  XM,YM, XT, YT,CRE,CRED, DIST, KILL 
END 


STOP 
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,«  V  *  -  V  *  ‘ 

%%  .  .V-V-V •/aV-*. 


n 


&££EHDIX  D 

PROGRAH  OF  THE  MANEUVERING  TARGET 


TITLE 

TITLE 

TITLE 

INTEG 

INTGER 

CONST 

CONST 

INITIAL 


BASIC  COMMAND  TO  LOS 
WITH  MANEUVERING  TGT 
***  YEUN,  J.Y.  *** 

RKSFX 

NPLOT,  KILL 
NPLOT  =  1  ,TCON=1.0 

VM=500. ,VT= 250, PI* 3.  1415  93,  KILL=0 


XTO  =  4000. 
YTO  =  1000. 
GAMT  =  PI 
F  =  0. 


***************************  ********************************* 
DERIVATIVE 

***************************  ****************************** *** 
NOSORT 

************************* 

***  TARGET  PARAMETERS  *** 

************************* 

VTX  =  VT*COS  (GAMT) 

VTY  =  VT*SIN]GAMT) 

XT  =  VTX*TIME  ♦  XTO 

YT  =  100*SIN(0.5*PI*TIME)  *YT0 

SIGT  *  ATAN2JYT.XT) 

IF  (TIME.  GE.TCON)  GO  TO  50 
************************* 

***  MISSLE  PARAMETERS  *** 

************************* 

*  PROGRAMMED  GUIDANCE  ************************************** 

* 

SIGM  *  SIGT  +0. 1 
VXM  =  VM  *  COS  (SIGM 
VYM  a  VM  *  SIN  (SIGM 
XM  =  I NTG  RL  (0 .  ,  VX  M 
INTGRL  (0.  .VYMi 
SQRT(XM**2  ♦  YM**2) 

=  SIGM  -  SIGT 
RM  *  SIN  (SIG  MT) 

DERI  V  (0.  ,CRE) 


YM 
RM  = 
SIGHT 
CHE  = 
CRED 


GO  TO  200 

***************************  ** 

*  ON-OFF  GUIDANCE  (BANG-BANG  CONTROL) 

******************  **—- -‘-‘- - — 

* 

CONTINUE 


t*******  ** 


********************* 
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G  =  150 
ACRED  a 
F  *  CRE 
U  =  -G 
CRED 


ABS  (CRED) 

♦  7cRED*AC RED) / (2*S) 
*  SIGN(1..P) 

INTGRL  (0 .0eU) 

D) 


^  0,U) 

CRE  =  INTGRL  (ORE, CRE 
ACRE  a  ABS  (CRE) 


RM  ■= 
A3  * 
SIGHT 
SIGM  « 


V  M*TI  ME 
CRE/RM 
»  ARSIN(A3) 

■  SIGT  ♦  SIGHT 


xm  =  rm  *  cos  (sigm) 

YM  =  R  M  *  SIN  SIGM) 

***************************  ******************  ************** 

200  CONTINUE 

* 

***************  ************  ******************************** 
***  MISSION  RESULT  ***  KILL  =  0  ;  IGT  MISSED 

**********************  KILL  =  1  ;  TGT  DESTROYED 

XD  1ST  =  XT-XM 
YDIST  =  YT-YM 

DIST  =  SQRT(XDIST**2  ♦  YDIST**2) 

IP  (DIST  .  LE.5)  KILL  =  1 

IF  (DIST  .GT.Sl  KILL  =  0 

IP  (XU  .GT.  (XT  +  30 )  )  CALL  ENDJOB 

SORT 

**************  ********************************************* 

********  OUTPUT  AND  PLOT  CONTROL  CARD  ********************* 
***************************  ******************************** 

SAMPLE 

CALL  DRWG  (1  ,  1,XM,  Y  M) 

CALL  DRWG  (1  ,2, XT,  YT) 

CALL  DRWG  (2, 1, TIME,  DIST) 

CALL  DRWG  3,1, TIME, CRE) 

CALL  DRWG  4,1, TIME, CRED) 

CALL  DRWG  (5,1,  CRE,  CRED) 

CALL  DRWG  (6 , 1,  TIME  ,  F) 

CALL  DRWG  (7, 1, TIME  ,U) 

TERMINAL 

CALL  ENDRW  (NPLOT) 

CONTRL  PINTIM=5.  9,  DELT=0.00  1,  DELS=0.  003 
PRINT  0.005,  XM,  YM,  XT,  YT,  CRE,  CRED,  DIST,  KILL 

END 
STOP 
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appendix  5 

PROGRAM  OP  THE  COMMAND  TO  PSEUDO-LOS 


TITLE  PSEDDO  -  LOS  COMMAND 

TITLE  WITH  IDEAL  RELAY 

TITLE  ***  YEUN,  J.Y.  *** 

INTEG  RKSFX 

INTGER  NPLOT  ,  KILL 

CONST  NPLOT=1  ,TCON=1.0 

CONST  VM  =500 . ,VT=250,PI=3.141593,  KILL=0 

INITIAL 

XTO  =  4000. 

YTO  =  1000. 

TG  =  0. 
f  =  0. 

GAMT  =  PI 

* 

***************  ************  **********************  *********** 
DERIVATIVE 

***************************  **********************  *********** 
NOSORT 

****  ****************  ***** 

***  TARGET  PARAMETERS  *** 

************************* 

VTX  =  VT*COS  (GAMT) 

VTY  =  V  T*SI  N  (GAMT) 

XT  =  VT  X*T  I  ME  ♦  XTO 
YT  =  VT Y*TI ME  +  YTO 
SIGT  =  ATAN2(YT.XT) 

XLOS  =  XT  +  VTX*TS 
YLOS  =  YT  +  VTY*T3 
SLOS  =  AT  AN  2  (YLOS,  XLOS) 

IP  (TIME.  GE  .  TCON)  GO  TO  50 
************************* 

***  MISSLE  PARAMETERS  *** 

************************* 

*  PROGRAMMED  GUIDANCE  ************************************** 

* 

SIGM  =  SLOS  +  0.1 
V MX  =  VM  *  COS  (SIGM) 

VMY  =  VM  *  SIN  (SIGM) 

XM  =  INTGRL  (0.  ,VM  X) 

YM  =  INTGRLfO.  .7MYJ 
RM  =  SQRT ( XM**2  +  YM**2) 

SIGMS  =  SIGM  -  SLOS 
CRE  =  RM  *  SIN  (SIGMS) 

CRED  =  DERIV  (0.  ,CR  E) 

GO  TO  2  00 

***************************  ** 

*  ON-OFF  GUIDANCE  (BANG-BANG  CONTROL)  ********************* 
****************** A ******** ** 

* 

50  CONTINUE 

* 

G  =  150. 

ACRED  =  A B S  (CRED) 

F  =  CRE  +  (CESD*ACRED)/ (2*G) 

U  =  -G  *  SIGN(1..P) 

CRED  =  INTG  RL  (0 . 5 ,  U) 

CRE  =  INTGRL  (CRE, CRED) 

* 

RM  =  V  M*TI ME 


A3  =  CRE/RM 
SIGHS  =  ARS IN ( A3) 

SIGH  =  SLOS  +  SIGHS 
VMX  =  VH  *  COS  (SLOS) 

VMY  =  VM  *  SIN  (SLOS 
XM  =  RM  *  COS  ( SI G  Ml 
Y M  =  RM  *  SIN  (SIGH1 

************************************************* *****1 

200  CONTINUE 

* 

**************  ******  ****  ***  ******************  ****  *****  3 

***  MISSION  RESO  LI  ***  KILL  =  0  ;  TGT  MISSED 

**********************  KILL  =  1  ;  TGT  DESTROYED 

XDIST  =  XT-XM 
YDIST  =  YT-YM 

DIST  =  SORT  (XDIST**2  +  YDIST**2) 

VC  =  SQRT ( (VTX-VMX) **2  + (VTY-VMY) **2) 

TG  =  DI  ST/VC 

IF  (DIST  .  L  E.  5)  KILL  =  1 

IF  (DIST  .  G  T.5)  KILL  =  0 

IF  (XM  .GT.  (XT  +  30)  )  l.ALL  ENDJOB 

SORT 

****************  ***  ********  ***************************< 
********  OUTPUT  AND  PLOT  CONTROL  CARD  ****************’ 

****  **********  *************  **********************  *****4 

SAMPLE 

CALL  DRWG  (1  ,  1,  XM,  Y  M) 

CALL  DRWG  (1  ,  2,  XT,  Y  T) 

CALL  DRWG  (2, 1,  TIME,  DIST) 

CALL  DRWG  (3,1, TIME, CRS) 

CALL  DRWG  (4  ,  1, TIME, CRED) 

CALL  DRWG  (5 , 1  ,CRE,  CRED) 

CALL  DRWG  (6 , 1,TIME  ,  F) 

CALL  DRWG  (7,1,  TIME  ,  U) 

TERMINAL 

CALL  ENORW  (NPLOT) 

CONTRL  FINTIM=5.  9,DELT=0.00  1  ,DELS=0.003 
PRINT  0.  1 ,T3, XM, YM  ,XT,YT, CRE , CRED , DIST, KILL 
END 


1PPJNDJX  P 

PROGRAB  OP  THE  BASIC  COMMAND  TO  LOS  WITH  TW0-L2VEL  RELAY 


TITLE 

TITLE 

TITLE 

INTEG 

INTGEH 

CONST 

CONST 

INITIAL 


COMMAND  TO  L .0. S 
TWO  LEVEL  RELAY 
YEUN,  J.Y.  ***** 


BASIC 
WITH 
***** 

RKSFX 
NPLOT, KILL 
NPL0T=1  ,TCO  N=1. 0 
VM*500.  ,VT=250,?I=3 


141593, KILL=0 


XTO  =  4  000, 
YTO  =  1000. 
GAMT  =  PI 
F  =0. 


***************************  ********************************* 
DERIVATIVE 

****  * *********  *****  ********  ***  ***************  ****  ********  *** 

NOSORT 

************************* 

***  TARGET  PARAMETERS  *** 

************************* 

VTX  =  VT*CO  S  (G  A  MT) 

VTY  =  VT*SI  N  ?GAMT) 

XT  =  VTX*TIME  +  XTO 
YT  =  VTY*TIME  ♦  YTO 
SIGT  =  AT  AN  2  (YT.  XT  ) 

IF  (TIME. GE. ICON)  GO  TO  50 
****************  ********* 

***  MISSLE  PARAMETERS  *** 

************************* 

*  PROGRAMMED  GUIDANCE  ************************************** 

* 

SIGM  =  SIGT  +  0.1 


VXM 
VYM  = 
XM  = 
YM  = 
RM  = 
SIGMT 
CRE  = 
CRED 


SIGM 


VM  *  COS 
VM  *  SIN, 

I NTGR  L  (0.  ,VXM) 
INTGRLfO.  .VYM 
SQRT(XM**2  ♦  YM**2) 
*  SIGM  -  SIGT 
RM  *  SIN  (SIGMT) 

=  DERI V  (0.  ,CR  E) 


GO  TO  2  00 
**** ************ *********** ** 

*  ON-OFF  GUIDANCE  (BANG-BANG  CONTROL) 


**** *******  ********** 


****************** ** ******* ** 
* 

50  CONTINUE 

* 


* 


G  =  150. 

ACRE  *  ABS  (CRE) 

ACRED  =  ABS  (CRED) 

IF  ( (ACRE*  ACRED)  .LT.  1.)  3  =  15 
F  *  CRE  ♦  JCRED*ACRED)  /  (2*3) 

U  =  -G  *  SIGN(1..F) 

CRED  =  INT G RL (0 . 0 ,  U) 

CRE  =  I  NTG  RL  (CRE,  C  RED) 

RM  =  V  M*TI  HE 
A3  =  CRE/RM 


SIGMT  =  ARSIN  (A3) 

SIGM  =  SIGT  ♦  SIGMT 
XM  =  RM  *  COS(SIGM) 

YM  =  RM  *  SIM  (SIGM) 

********************* ****** ********************** ********** 

200  CONTINUE 

* 

***  ***********  *************  ******************  ****  ********** 
***  MISSION  RESULT  ***  KILL  =  0  ;  TGI  MISSED 

**********************  KILL  =  1  ;  TGT  DESTROYED 

XDIST  =  XT-XM 
YDIST  =  YT-YM 

DIST  =  SQR  T (XDIST  **2  +  YDIST**2) 

IP  (DIST  .  LE.5)  KILL  =  1 

IF  (DIST  .GT.5i  KILL  =  0 

IF  (XM  .3T.  (Xx+30  )  )  CALL  ENDJOB 

SORT 

*************************************************  ********** 

********  OUTPUT  AND  ELOT  CONTROL  CARD  ********************* 
***************************  ******************************** 


SAMPLE 

CALL  DRWG  (1 
CALL  DRWG  1 
CALL  DR  WG  (  2 
CALL  DRWG (3 
CALL  DRWG  i  4 
CALL  DRWG  i  5 
CALL  DRWG  I  6 
CALL  DRWG  (7 

TERMINAL 


CALL  ENDRW  ( NPLOT) 

CONTRL  FINTIM® 5 .  9,DSLT=0.00  1,DELS=0.  003 
PRINT  0.005#XM,YM,XT,YT,CRE,CRED,DIST,  KILL 
END 


STOP 
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APPENDIX  G 

PHOGHAH  OF  THE  MANEUVERING  TARGET  0ITH  TWO-LEVEL  RELAY 


TITLE  BASIC  COMMAND  TO  LOS 

TITLE  WITH  MANEUVERING  IGT 

TITLE  ***  YE  UN  ,  J.Y.  *** 

INTEG  RKSFX 

INTGER  NPLOT  ,  KILL 

CONST  NPLOT=1  ,TCON=1.0 

CONST  VS  =500.  ,VT=  250,  PI  =3.141593,  KILL=0 

INITIAL 

XTO  =  4000. 

YTO  =  1000. 

GAMT  =  PI 
F  =  0. 

♦ 

***************************  ********************************* 
DERIVATIVE 

***************************  ********************************* 

NOSORT 

************************* 

***  TARGET  PARAMETERS  *** 

************************* 

VTX  =  VT*CO  S  (GAMT) 

VTY  =  VT*SIN?GAMT 
XT  =  VT  X*TI  ME  +  XTO 
YT  =  1 0 0*SI N  (0. 5*P I*TIME) +YT0 
SIGT  =  AT AN2 JYT.XT) 

IF  (TIME.GE  .ICON)  GO  TO  50 
****************  ********* 

***  MISSLE  PARAMETERS  *** 

************************* 

*  PROGRAMMED  GUIDANCE  ************************************** 

* 


SIGM  =  SIGT+0.1 
VXM  =  VM  *  COS(SIGM) 
VYM  =  VM  *  SIN  (SIG  M 
XM  =  INTGRL  (0.  , VXM 
YM  =  I  NTGRL  (0.  ,  VI  M 
RM  =  SQRT  ( XM**2  «■  YD! 
SI  GMT  =  SIGM  -  SIGT 
CRE  =  RM  *  SIN  (SIS  MT) 
CRED  =  DERI V  (0 .  ,CR  E) 


M**2) 


GO  TO  2  00 

***************************  ** 

*  ON-OFF  GUIDANCE  (BANG-BANG  CONTROL) 
***************************  ** 

* 

50  CONTINUE 


********************* 


G  =  150. 

ACRE  *  ABS (CRE) 
ACRED  *  ABS  (CRx 
IF  (  (ACR  E+ ACRED) 
F  =  CRE  ♦  (CREI 
U  =  -G  *  SIGNM 
CRED  =  INTGRL  (C 
CRE  =  INTGRL  (Cl 
ACRE  =  ABS  (CRE) 


CRE) 

(CRED) 

RED)  .LT .  1.)  G  =  15. 
CRED*ACRED)  /  (2*G) 


CRE, CRED) 


RM  =  V  M*TI  ME 
A3  =  CRE/RM 
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SIGHT  =  4RSIN  (A3) 

SIGM  =  SIGT  ♦  SIGMT 
IN  =  RM  *  COS  (SIGM) 

YM  =  RM  *  SIN  (SIGH) 

***************************  ******************************** 

200  CONTINUE 

* 

***************************  ******************************** 
***  MISSION  RESULT  ***  KILL  =  0  ;  TGT  HISSED 

**********************  KILL  =  1  ;  TGT  DESTROYED 

XDIST  «  XT-XM 
YDIST  *  YT-YM 

DIST  =  SQRT  <XDIST**2  ♦  YDIST**2) 

IF  (DIST  .LE.5)  KILL  =  1 

IF  (DIST  .  GT.Sj  KILL  =  0 

IF  (XM  .GT.  <XT+3Q))CALL  ENDJOB 


SORT 

***************************  ******************************** 

********  OUTPUT  AND  PLOT  CONTROL  CARD  ********************* 
***************  ************  ******************************** 

SAMPLE 

DR  WG  (1  ,  1 ,  XM.  X  M) 

DRWG  1 , 2,  XT,  Y  T) 

DR  WG  i  2 , 1, TIME  , DIST) 

DRWG  3 , 1, TIME, CHE) 

DRUG  I  4 ,1,  TIME  .CRED) 

DRWG  ,  5,  1 ,  CRE,  CRED) 

DR  WG  (6,1,  TIME  ,  F) 


TERMINAL 


CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 


CONTRL 

PRINT 

END 

STOP 


DRWG  (7;i;timeJu{ 
NPLOT) 


CALL  ENDRW  (1 

FINTIM*  5.  9,  DE  LT*0 .00  1  ,  DBLSs0.  0  03 
0. 005, XM, YM, XT, YT, CRE, CRED, DIST, KILL 
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APPENDIX  B 

PBOGB&H  OP  THE  COMMAND  TO  PSEUDO-LOS  HITH  TNO-LEVEL  BELAY 


TITLE  PSEODO  -  LOS  COMMAND 

TITLE  ilTH  THO-LEVEL  BELAY 

TITLE  *****  YEUN,  J.Y.  ***** 

INTEG  RKSFX 

INTGEB  NPLOT ,  KILL  ^  A 

CONST  NPLOT=1  .TCON*1 .0 

CONST  VM-500.  JVT=250,  PI *3.  141593,  KILL=0 

INITIAL 

XTO  =  4000. 

YTO  =  1000. 

TG  *  0. 

GAMT  =  PI 
P  =  0. 

***************************  ********************************* 

***************************  ********************************* 

NOSOHT 

************************* 

***  TARGET  PARAMETERS  *** 

************************* 

VTX  *  VT*COS  (GAMT) 

VTY  *  V T*SI  N JG AMlj 
XT  =  VT X*TI ME  ♦  XTO 
YT  *  VT Y*TIME  ♦  YTO 
SIGT  =  AIAN2  (YT.XT) 

XLOS  *  XT  ♦  VTX*TG 
YLOS  »  YT  ♦  VTY*TG 
SLOS  *  ATAN2JYLOS,  XLOS)  „ 

IF  (TIME.  GE.TCON)  GO  TO  50 
************************* 

***  MISSLE  FAB AMETEBS  *** 

************************* 

*  PROGRAMMED  GUIDANCE  ************************************** 

* 

SIGM  a  SLOS  ♦  0.1 
VHX  =  VM  *  COS  (SIGM) 


RM  »  SQRT(XM**2  ♦  YM**2) 
SIGHS  »  SIGM  -  SLOS 
CRE  *  RM  *  SIN  (SIG  NS) 

CRED  *  DERIV(0.  ,CRE) 

GO  TO  200 

***************************  ** 

*  ON-OFF  GUIDANCE  (BANG-BANG  CONTROL) 
***************************  ** 

* 

50  CONTINUE 


********************* 


G  =  150. 

ACRE  *  ABS(CBE) 
ACRED  *  ABSjCRI 


IP  ((ACRE*  ACRED)  .LT.  1.)  G  * 
P  »  CRE  *  (CREu*AC RED) / (2*G) 
U  *  -G  *  SIGN(1.,P) 


»  15. 


U  »  -G  *  SIGN 
CRED  »  INTGRL 
CRE  *  INTGRL ( 


jo  1 6  fin 

CRE,l REI 


RM  =  V  M*TI  ME 
A3  =  CRE/RM 
SIGMS  »  ARSIN (A3) 

SIGH  *  SLOS  ♦  SIGMS 
VMX  =  VM  *  COS(SLOS) 

VMY  =  VM  *  SIN  (SLO S 
XM  =  R  M  *  COStSIGM 
Y M  =  HM  *  SIN  (SIG  M 

****  *** **** *****  **  ****  ******************************** 

200  CONTINUE 

* 

********************************************************* ** 

***  MISSION  RESULT  ***  KILL  =  0  ;  TGT  MISSED 

**********************  KILL  =  1  :  TGT  DESTROYED 

XDIST  =  XT-XM 
YDIST  =  YT-YM 

DIST  =  SQRTjXDIST**2  +  YDIST**2) 

VC  =  SORT ( (VTX-VMX) **2+ (VTY-VMY) **2) 

TG  *  DI  ST/VC 

IF  (DIST  .LE.5)  KILL  =  1 

IF  (DIST  .GT.5j  KILL  =  0 

IF  (XM  .GT.  (XT+30)  )  CALL  ENDJOB 

SORT 

***************************  ******************************** 

********  OUTPUT  AND  PLOT  CONTROL  CARD  ********************* 
********************************************************* ** 


SAMPLE 


TERMINAL 


CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 


DRWG 

DRWG 

DRWG 

DRWG 

DRWG 

DRWG 

DRWG 

DRWG 


(1  ,  1r  XM,  Y  M) 

ri  #  2f  XT#  yt{ 

2  f  If  TIME  r  DI 
3 , If  TIME , CR 
4.1. TIME. CR 


)IST) 


CALL  ENDRW  (NPLOT1 

CONTRL  FINTIM*5.  7,  DELT=0.D0  1  .DELS=0.003 
PRINT  0. 1,TG, XM, YM  fXT,YT, CRE ,CRED, DIST , KILL 
END 


STOP 
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APPENDIX  I 

PROG BAH  OF  THE  BASIC  COHHAND  TO  LOS  KITH  SATURATION  CONTROL 


TITLE  BASIC  COHHAND  TO  L.O.S 

TITLE  MISSLE  CONTROL 

TITLE  WITH  SATURATION  CONTROL 

INTEG  RKSFX 

INTGEB  NPLOI,  KILL.  CUR 

CONST  NPLOT=1 ,TCO N=1 . 0, C UR= 1 

CONST  VH=500 . ,V T= 250, PI  = 3 . 1 41 5 93,  KILL=0 

PARAM  H  =  10. 

INITIAL 

XTO  =  4000. 

YT0  =  1000. 

GAHT  »  PI 

F  =  0. 

* 

***************************  **********************  *********** 
DERIVATIVE 

***************************  ********************************* 
NOSORT 

************************* 

***  TARGET  PARAHETERS  *** 

****  ********************* 

VTX  a  VT*CO S  (GAHT) 

VTT  »  VT*SIN?GAMTi 
XT  *  VTX*TIHE  ♦  XTO 
YT  =  VT Y*TIHE  ♦  YT0 
SIGT  =  AIAN2  (YT.XT) 

IF  (TIHE.GS.ICOlf)  GO  TO  50 
****************  ********* 

***  HISSLE  PARAHETERS  *** 

************************* 

*  PROGRAMHED  GUIDANCE  ************************************** 

* 

SIGH  =  SIGT  +  0.1 
VXH  a  VH  *  COS  (SIGH) 

VYM  =  VH  *  SIN  (SIG  H 
XM  a  INTGRL  (0.  ,  VX  H 
YH  *  INTGRL  (0.  .  VY  H 
RH  =  SQRT(XH**5  ♦  IH**2) 

SIGHT  =  SIGH  -  SIGT 
CHE  =  RH  *  SIN  (SIGHT) 

CRED  =  DERI V  (0.  ,CR  E) 

GO  TO  2  00 

******************** ******* ** 

*  ON-OFF  GUIDANCE  (BANG-BANG  CONTROL)  ********************* 
***************************** 

* 

50  CONTINUE 


H**2) 


CONTROL)  ********************* 


G  a  150. 

ACRED  a  ABS(CRED) 

F  a  CRE  *  ( CRED* ACRED) 
U  =  -G  *  SI GN ( 1 . # P ) 

IF  (ABS(F)  .LT.  H)  U= 
CRED  »  INTG  RL (0. 0 .  U) 


S  {FI  .LT.  H)  U= 
INTGRLfO.O,  U) 
INTGR  L  (CRE,  C  RED) 


/ (2*G) 
-G*F/H 


RH  a 
A3  * 
SIGHT 


V  H*TIHE 
CHE/RH 
a  ARSIN 


SIN  (A3) 
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SIGH  =  SIGT  ♦  SIGMT 
XM  =  R  M  *  COS(SIGM) 

YM  =  RM  *  SIN (SIG  Mj 

*************************** ****************** *** *  ********** 

200  CONTINUE 

* 

***************************  ******************  ************** 


***  MISSION  RESULT  *** 

KILL  = 

0 

;  TGI  MISSED 

********************** 

KILL  = 

;  TGT  DESTROYED 

XDIST  =  XT-XM 

YDISI  =  YT-YM 

DIST  =  SORT  (XDIST  **2 

* 

YDIST**2) 

IF  (DIST  .  LE.5) 

KILL 

s 

1 

IF  (DIST  .  GT.5) 

KILL 

= 

0 

SORT 

*********************************************************** 

********  OUTPUT  AND  PLOT  CONTROL  CARD  ********************* 
*********************************************************** 

SAMPLE 

CALL  DRUG  (1  ,CUR,XM,YM) 

CALL  DRHG  (2  ,CUR  ,TI  ME,  DIST) 

CALL  DRWG  (3  ,CUR,TIME,  CRE) 

CALL  DRWG  (4  ,CUR,CRE.CRED) 

CALL  DRWG  (5  ,CUR,riME,  F) 

CALL  DRWG  (6  ,CUR,TIME,  U) 

TERMINAL 

IF  (CUR  .EQ.  3)  CALL  ENDRW  (NPLOT) 

CU  R  *  CUR  +  1 

CONTRL  FINTIM=  5. 65 , D ELT=0. 0 0 1 , DELS=0 . 00 3 
PRINT  0.005,XM,YM,  XT,  YT,  CEE,  CRED,  DIST,  KILL 
END 

PARAM  M  =  5. 

END 

PARAM  M  =  1. 

END 

STOP 
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APPENDIX  J 

PROGBAB  OF  THE  MANEUVERING  TARGET  BITH  SATURATION  CONTROL 


TITLE  BASIC  COMMAND  TO  LOS 

TITLE  (MANEDVERING  TGTf 

TITLE  WITH  SATURATION  CONTROL 

INTEG  RKSFX 

INTGER  NPLOT,  KILL, CUR 

CONST  NPIOT=1  #TCON=1.0,CUR=1 

CONST  VM=500  .  ,VT=250,PI  =  3.1415  93,  KILL=0 

PARAM  M  =  10. 

INITIAL 

XTO  =  4000. 

YTO  =  1000. 

GAMT  =  PI 
F  *  0. 

* 

***************************  **********************  *********** 
DERIVATIVE 

***************************  ********************************* 

NOSORT 

************************* 


***  TARGET  PARAMETERS  *** 

************************* 

VTX  ■  VT*COS  (GAMT) 

VTY  =  VT*SIN]GAMTj_ 

XT  »  VT X*TIME  ♦  XTO 

YT  =  100*SIN(0.5*PI*TIME) +YT0 

SIGT  =  ATAN  2jYT,XT) 

IF  (TIME.  GE. ICON)  GO  TO  50 
************************* 

***  MISSLE  PARAMETERS  *** 

************************* 

*  PROGRAMMED  GUIDANCE  ************************************** 

* 


SIGH  a  SIGT  +0. 1 
VXM  a  VM  *  COS  (SIGH 
VYM  =  VM  *  SIN  (SIS  M 
XM  a  I  NT G R L  (0  .  *  V X  M 
YM  *  INTGRL(0.  .VYM 
RM  »  S  QRT (XM**2  ♦ 
SIGHT 
CRE  = 


JL  HA  VJU  IV  «  |  f  A  lil 

SQRT(XM**2  ♦  YM**2) 
a  SIGH  -  SIS  T 


CRE  =  RM  *  SIN  (SIGHT) 
CRED  *  DERIV(0.  ,CRE) 
GO  TO  2  00 

***************************  ** 


*  ON-OFF  GUIDANCE  (BANG- BANG  CONTROL) 

***************************  ** 


********************* 


CONTINUE 


G  =  150. 

ACRED  a  ABS(CRED) 

F  »  CRE  ♦  (CRED*ACRED)/ (2*3) 
U  =  -G  *  S IGN  (1mF) 

IF  (ABS  (F)  .  LT.M)  Ua-G*F/M 
CRED  »  INTGRL  (0.0  ,U) 

CRE  a  INTGRL(CRE,  CRED) 


ACRE  *  ABS(CRE) 


RM  = 
A3  = 


V  M*TI  ME 
CRE/RH 
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SIGHT  =  ARSIN  (A3) 
SIGH  =  SIGT  ♦  SIGHT 


XM  = 
YM  = 


[SIG  Ml 


*********************************************************** 

200  CONTINUE 

* 

***************************  ******************************** 


***  MISSION  RESULT  ***  KILL  = 
**********************  RILL  = 
XDIST  *  XT-XM 
YDIST  =  YT-YM 
DIST  =  SQR  T (XDI5T  **2 


MISSED 

DESTROYED 


LE.5) 
i  GT. 5 


YDIST**2) 

0 


IF  (DIST  .LE.5)  KILL  =  1 
IF  (DIST  .  GT. 5)  KILL  =  0 

* 

SORT 

***************************  ******************************** 

********  OOTPUT  AND  PLOT  CONTROL  CARD  ********************* 
***************************  ******************************** 

SAMPLE 

CALL  DRHG  (1  ,CUR,XM,YM) 

CALL  DRHG  i  2  ,CUR,TI  ME, DIST) 

CALL  DRHG  3  ,CUR,TI  ME,  CRE) 

CALL  DRHG  (4  ,CUR,CRE,CRED) 

CALL  DRHGi5,CUH,riME,  F) 

CALL  DRHG  (6  ,CUR,TI  ME,  U) 


TERMINAL 


CALL 

CALL 

CALL 

CALL 

CALL 

CALL 


DRHG 

DRHG 

DRHG 

DRHG 

DRHG 

DRHG 


IF  (CUR  .  EQ  .  3)  CALL  ENDRH  (NPLOT) 

COR  =  CUR  ♦  1 

FINTIM=5.  65,DELT=0.001,DELS=0. 003 
0.005,XM,YM, XT, YT,  CRE, CRED, DIST, KILL 

M  =  5. 

M  *  1. 


CONTRL 

PRINT 

END 

PAR  AM 

END 

PARAM 

END 

STOP 


APPENDIX  K 


PROG RAH 

TITLE 

TITLE 

TITLE 

INTEG 

INTGEH 

CONST 

tm  '  ,  ’ 

CONST 

JS,/’ 

PARAM 

INITIAL 

TITLE  PSEDDO  -  LOS  COMMAND 

TITLE  WITH  SATURATION  CONTROL 

TITLE  *****  YEUN,  J.Y.  ***** 

INTEG  BKSFX 

INTGEH  NPLOT, KILL, CUR 

gSI?  ?Iii5s:;5?^lJ;?i£03!TJ,593.KiLL=o 

PARAM  M  =  10. 

INITIAL 

XTO  *  4  000. 

YTO  *  1000. 

TG  =  0. 

GAMT  *  PI 
F  =0. 

***************************  ********* 1  ***********  *********** 
SI************************* ********************************* 

NOSORT 

************************* 

***  TARGET  PARAMETERS  *** 

************************* 

VTX  *  YT*CO  S  (GAMT) 

VTY  =  V  T*SI  NjfGAMT) 

XT  =  YT X*TI ME  ♦  XTO 
YT  =»  VTY*TIME  ♦  YTO 
SIGT  =  ArAN2jfrr,XT) 

XLOS  *  XT  ♦  vTX*T3 
YLOS  =  YT  ♦  VTY*TG  , 

SLOS  =  ATAN2(YLOS,XLOS) 

IF  (TIME.  GE.TCON)  GO  TO  50 
************************* 

***  MISSLE  PARAMETERS  *** 

**PROGRAMHED*GUIDANCE******  ********************************* 
* 


SIGH  =  SLOS 


YHX  * 
YHY  = 
XH  * 
YM  * 
RM  » 
SIGHS 
CRE  = 
CRED  = 
GO  TO 


YM  *  COS  (SIGM 
VM  *  SIN  (SI 3  M) 

I  NT  G  R  L  (0 .  ,  7  M  X 
INTGBL  (0.  ,7M  YV 
SQRT  (XM**2  ♦  YK 
*  SIGM  ~  SLOS 
RM  *  SIN  (SIGHS) 
DERI  7  (0.  ,CR  E) 
200 


M  **  2) 


***************************  ** 

*  ON-OFF  GUIDANCE  (BANG-BANG  CONTROL) 
***************************  ** 

* 

50  CONTINUE 
*  G  =  150. 

)C2E?8I  }M*fflacUDl/(2.3| 
ip'disirf  P|  -  -s  *  p 


********************* 


KVJ 


RM  =  7  M*T I ME 
A3  =  CRE/RM 
SIGMS  =  A  RS  IN  ( A3) 

SIGH  =  SLOS  ♦  SIGN  S 
VMX  =  VM  *  COS  (SLOS 
VMY  =  VM  *  SINiSLOS 
XM  =  RM  *  COS  (SIGM 
YM  =  RM  *  SIN  i' SIS  M 

***************************  **********************  ******** ** 
200  CONTINUE 

* 

***************************  ******************************** 
***  MISSION  RESULT  ***  KILL  =  0  ;  TGT  MISSED 

**********************  KILL  =  1  ;  TGT  DESTROYED 

XDIST  =  XT-XM 
YDIST  =  YT-YM 

DIST  =  S  QRT  (XDIS T  **2  ♦  YDIST**2) 

VC  =  SORT  (  (  VTX-VMX)  **2+  (VTY-VMY)  **2) 

TG  =  DIST/VC 

IP  (DIST  •  LE.5)  KILL  =  1 

IF  (DIST  .GT.5)  KILL  =  0 

SORT 

*********************************************************** 

********  OUTPUT  AND  PLOT  CONTROL  CARD  ********************* 
****  ***********  ****************************  **  ****  ********  ** 

SAMPLE 

CALL  DRWG (1 ,CUR,XM ,YM) 

CALL  DRUG  I  2  ,CUR,TI  ME, DIST) 

CALL  DRWG  (3  ,  CUR,TI  ME,  CRE) 

CALL  DRWG  (4  ,CUR,CRE,CRED) 

CALL  DRWG  i  5, CUR, TIME,  F) 

CALL  DRWG  I  6  ,CUR,TI  ME,  U) 

TERMINAL 

IF  (CUR  .EQ.3)  CALL  ENDRW  (NPLOT) 

COR  *  COR  +  1 

CONTRL  FIN  TIM*  5.  65  ,  DELT=0.  0  01  ,DELS*0.003 
PRINT  0.0C5,TG,XM, YM, XT, Y T , CRE, CRED, DIST , KILL 

PARAM  M  =  5. 

END 

PARAM  H  *  1. 

END 

STOP 
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